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This  thesis  describes  the  development  of  a  new  method 
for  controlling  the  error  in  model  following  control  syst 
The  treatment  is  for  first  order,  linear  or  nonlinear, 
time  varying  or  time  invariant  systems  with  additive 
(linear)  control.   The  errors  controlled  are  assumed  to 
have  arisen  from  external  disturbances  or  from  differences 
in  the  initial  conditions  of  the  plant  and  the  model. 

The  theory  introduced  is  called  constant  error  model 
following  control.   This  paper  describes  the  theory  as  an 
outgrowth  of  attempts  to  control  a  plant  by  feedback  of  an 
error  between  the  plant  and  model  with  the  error  specified 
to  be  constant.   From  this,  it  is  shown  that  a  model  may  be 
followed  with  arbitrary  error.   The  central  result  is 
that,  given  some  error,  one  can  find  another  model  (control 
model)  which,  if  followed  with  this  arbitrary  error,  will 
guide  the  plant  state  trajectory  back  to  that  of  the  model. 
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I .   INTRODUCTION 

The  theories  presented  in  this  thesis  were  developed 
in  the  course  of  researching  control  techniques  for  use  in 
aircraft  departure  prevention.   The  aim  was  to  apply  model 
following  control  methods  to  the  nonlinear,  time  varying 
and  coupled  dynamic  equations  of  motion  of  an  aircraft  at 
critical  combinations  of  angle  of  attack  and  angular  rates. 
When  it  was  seen  that  a  new  development  in  model  following 
control  had  evolved  from  the  research,  it  became  the  sole 
subject  of  the  thesis. 

Model  following  control  is  concerned  with  techniques 
which  cause  a  physical  plant  to  behave  as  much  like  a  model 
as  possible.   Motyka  [Ref.  1]  accomplishes  this  by  solving 
the  plant  state  equations  for  their  controls  and  then 
substituting  expressions  for  the  model  states  and  state 
rates  into  these  equations.   That  is,  assuming  that  the 
plant  is  defined  by  the  linear,  small  perturbation  constant 
coefficient  differential  equation: 

X   =  f-,  X   +  g,   u  (1) 

p    1  p   °1  p 


where   x~  is  the  state  of  the  plant 
P  F 

u    is  the  control 

P 


and  that  a  model  is  given  by  the  corresponding  differential 
equation : 

x   =  f0  x   +  o      u  (21 

m    2   m   62   m  ^    J 

It  is  desired  that   x   =  x    and   x   =  x 

m    p        m    p 

The  plant  control  which  achieves  this  can  be  determined  by 
substituting  the  desired  relationships  into  the  plant 
equation : 


x   =  f,  x   +  g,  u  (3) 

m    1   m   tol   p  ^  J 


And  solving  for   u   .   The  result  is  an  expression  for  the 
plant  inputs  which  make  the  plant  state  equal  to  that  of 
the  model. 

Problems  with  this  method  arise  when  equations  (1) 
and  (2)  are  vector  equations  with  fewer  controls  than  states 
to  be  controlled.   Also,  if  equation  (1),  which  is  a 
mathematical  description  of  a  real,  physical  plant,  fails 
to  describe  that  plant  accurately,  errors  may  be  introduced 
into  the  response  of  the  system.   Finally,  errors  will 
occur  if  the  plant  and  the  model  do  not  have  the  same 
initial  conditions,  or  equivalen tly ,  in  the  presence  of 
external  disturbances. 

The  first  two  of  these  three  problems  are  treated  by 
Moytka  [Ref.  1]  and  Rynaski  [Refs.  2  and  3].   They  are 
beyond  the  scope  of  this  paper.   The  theories  presented 


herein  apply  specifically  to  the  problem  of  reducing  the 
system  error  in  the  third  case.   The  method,  called  constant 
error  model  following  control,  is  developed  and  demonstrated 
for  nonlinear,  time  varying,  first  order,  single  input - 
single  output  systems  with  additive  (linear)  control. 
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II.   PROBLEM  DESCRIPTION 

For  the  purposes  of  this  discussion,  it  is  assumed 
that  a  plant  (physical  process)  is  given,  and  that  it  is 
completely  described  by  a  first  order  differential  equation 
with  additive  control  of  the  form 


x„  =  f,  (x„  ,  t)  +  a      u  (1) 

p     1  v  p  '  J         &1   p  y    J 


where   x    is  the  state  of  the  plant 
P  v 

t    is  the  variable  time 

u    is  the  control  input  to  the  plant 
p  v  v 

f,  is  (in  general)  time  varying  and  nonlinear  in  x 

g,  is  constant 

J-  • 

It  is  desired  that  the  response  of  the  plant  be  modified  in 
some  way.   This  modified  response  is  completely  described 
by  the  mathematical  model 


x   =  £-  (x   ,  t)  +  g9  u  (2) 

m     2  ^  m        &2   s 


where   x    is  the  state  of  the  model  which  has  the  desired 
m 

response 

u    is  the  control  input  to  the  system  which  incor- 
porates the  model  and  the  plant,  and 

f9   is  (in  general)  time  varying  and  nonlinear  in 
x   ,  and   f9  f    f- 


cr 


2   is  constant,  g9  f    g-. 
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III.   SYSTEM  DESCRIPTION 


A  parallel  model  following  system  is  chosen 


U- 


/A0D6L. 


X 


M 


U£_ 


►  rl^NT 


Ve. 


U* 

.Ka 

FIGURE  1.   PARALLEL  MODEL  FOLLOWING  CONTROL  SYSTEM 
where   e   is  the  magnitude  of  the  system  error: 

e  =  I  x   -  x 

1  m    p ' 

K    is  the  adaptive  gain  applied  to   e  to  generate 

uT, 
uv   is  the  (additive)  modification  to  the  system 


lK 


control ,  u 


In  this  system,  the  steady-state  error  (e   )  can  never  be 
zero  unless  the  plant  and  model  have  identical  responses 
to  the  same  input,  since 

e    =0  implies  that  uv   =    0  implies  that   u   =  u 


ss 


P     s 


and 


e    =0  implies  that   x_„  =  x 
ss       *  mss     pss 
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Since  this  condition  serves  no  useful  purpose  unless  the 
parameters  of  the  plant  or  model  are  adjustable,  we  intro- 
duce the  idea  of  a  constant  (non-zero)  error. 
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IV.   CONSTANT  ERROR  MODEL  FOLLOWING  CONTROL 

A.   CONSTANT  ERROR 

Without  changing  the  system,  the  error  (e)  is  specified 
to  be  fixed  at  some  arbitrary  value  (e): 


U* 


$H  MooeL 


*M 


l~>&: 


uj 


PuAMT 


V? 


K* 


4.  e«£ 


FIGURE  2.   CONSTANT  ERROR  MODEL  FOLLOWING  CONTROL  SYSTEM 

We  wish  to  determine  the  gain,   K*  ,  which  will  insure  that 
this  condition,  once  established,  will  be  maintained  as 
shown  in  Figure  3. 


FIGURE  3.   PLANT  AND  MODEL  STATE  TRAJECTORIES 

The  state  trajectory  of  the  plant  follows  that  of  the  model 
with  constant  difference   e  . 
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We  have  defined 

e  =  |x   -  x  |  (4) 

1  m    p1  v  } 

A 

so  that   £  =  |x   -  x  *  I  =  constant, 

where   x  *  is  the  plant  output  for   e  =  z    .      Similarly, 
P 

denote  the  value  of   K  which  maintains  the  error   constant 
at   e  =  £   as   K*   and  the  resulting  control  as   u*  . 
W e  now  have 


e  =  |x   -x*|=(x   -x*)sgn(x   -x*) 
1  m    p1    ^  m    py&vm    p 


x   sgn  fx   -  x  *)  -  e 
*    m   *   ^  m    p 

SO  X   *   =  t : t-ar^ 

p       sgn  (x   -  x  *J 


or 


x  *  =  x   -£sgn(x   -  x  *)  (5) 

p     m      &    m    p 


and  z    =  0  =  (xm  -  xp*)  sgn  (xm  -  xp*) 


r°_    [scrn    fx      -    x    *)]    =    0      since      (x      =   x   *)       is    assumed 

Mt    L    °      *■    m  P  m  P 


constant  ) 


so  x    *    =    x 

D  m 
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The  result  is  that,  if  the  plant  state  is  exactly  +_  e 

from  the  model  state,  we  may  express   x  *   and   x  *   as 

P         P 

functions  of   x   ,  x   and  +  e  .   The  values  of   x   ,  x 

mm—  mm 

and   e   are  known  for  any  time  and  system  control  because 
they  have  been  specified. 

Now  note  from  the  system  description  (Figure  2)  that 


u  *  =  u   +  u  * 
p     s     K 


or  u  *  =  u  *  -  u  (7) 


and  u  *  =  K  *  e  (8) 

K     a  *•  J 


Equations  (1)  and  (2)  may  be  solved  for   u  *   and   u,  : 


1 


x   -  (f,  (x  *  ,  t) 

u  *  =  _J2 1    P  (9) 

P  gx 


x   -  £-  (x   ,  t) 

=  _m 2 m (1Q) 

s  go 


By  substituting  equations  (5)  and  (6)  into  equation  (9) 


x   -  f,  [(x   -  e  sgn  (x   -  x  *))  ,  t] 

u  *  =  m  I  m I    m    P     -      rii) 

P  gi 


If   U   and   u    are  not  additive  controls  as  assume 
P        s 

and  if  the  state  equations  may  be  solved  for   u    and   u   , 
the  results  which  follow  are  still  valid.       H 
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Here   we    define 


V  (xm  -  e  •  l  )  k  fi  tt-v  -  e  sEn  c*m  -  v»  >  '] 


m         v  m    p 


(12) 


■  fi  (XP*  >  v 


x   -  f ,  *  (x   ,  £  ,  t) 

up*  -  -£ L_liU —         (13) 


From  equations  (7),  (8),  (10)  and  (13) 
K*  e  =  u*  =  u*   -  u 


X       f  *   (X  ,   £   ,   t)      X    "   f0   (X    ,   t) 

K*  E  =   m  -  1 m  =   B     2    5  r141 

°1  &2 

x   -  f  *  (x  ,  £  ,  t)    x   -  f0  (x   ,  t  ) 

K*    _  _m 1   2  m  .  _JB 2  ^  m  /--.r-j 

°1  &2 


Equation  (15)  gives  us  the  gain  which  will  insure  that, 

if   x   =  x   +  £   initially  (Figure  4) ,  the  proper  control 
p     m  -  /  v.   6      j  ,  f      f 

will  be  applied  to  maintain  it  there  if  the  system  is 
undisturbed.   It  is  expressed  solely  in  terms  of  the  model 
state  and  state  rate,  the  specified  error,  and  the  func- 
tional relationships  which  define  the  plant  and  model 
responses . 
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FIGURE  4.   PLANT  AND  MODEL  STATE  TRAJECTORIES 

The  use  of   K*   gain  is  illustrated  in  the  following 
example : 

Example  1:   Linear,  first  order,  time  invariant  model  and 
plant . 


x   =  a  x   +  u 

P      P    P 


x   =  b  x   +  u 

m      m    s 


a  f   b 


we  have 


x*  =  ax*  +  u* 

P       P     P 
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then 


u*    =   x*-ax* 

P  P  P 


u        =    x      -    b   x 
s  m  m 


u   *    -    u 


K    *    =   _E 1 

a  z 


(x      -ax*)-fx      -bx) 
p  p    *         K   m  nr 


From   equations    (5)    and    (6) 


x 

K   *    =   -2L 
a 


a[x      -    e    sgn    fx      -    x   *) 1    -    x      +   b   x 

L   m °  m p      J m m 


or 


(b    -    a)    x      +   a  £    sen    (x      -    x   *) 

K   *    = m m  p 

a  e 


(16) 


The  system  is  as  shown  in  Figure  5 


UL 


Xm 


uu" 


&-£•■* 


«£{ 


*£ 


CL 


FIGURE  5.   CONTROL  SCHEMATIC  FOR  EXAMPLE  1 
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The  system  was  simulated  on  the  IBM  3033  using  the 
Continuous  System  Modeling  Program  (CSMP) .   The  following 
values  were  used: 

a  =  +  0.5        (unstable  impulse  response) 

b    =    -    1.0 

e    =    0.05    2 

x    (0)    =    0 

mv    -* 


x    (0)    =    0.1 
P 


XDCT  (PLPNT!  =  0.5*X  (PLANT)  +U  (PLANT) 
B  XDGT(MQCEL)=-:.*X  (MODEL)  +  U  (SYSTEM) 
K  PLANT     I. 

EPSILQN 


LEGEND 
XM  □ 

XP  O 


^.00 


0.5O  1.00 


I. SO  2.00  2. SO  3. OO  3.50  4.00  U.50  5.00 

TIME 


FIGURE  5a.   SIMULATION  RESULTS  FOR  EXAMPLE  1 


This  value  of   e   was  chosen  so  that  the  error 
following  could  be  shown  graphically.   The  system  performs 
similarly  for  any  value  of   e  >  0  . 
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The  input  (u  )  was  a  unit  step  at   t  =  0  .   A  tin 


me 


history  of  the  system  response  for  the  first  five  seconds  is 
shown  in  Figure  5a. 

B.   STABILITY  OF  THE  SYSTEM  USING   K*e   CONTROL 

The  question  to  be  answered  is:   Will  the  gain   K* 

satisfactorily  control  the  plant  if   e  f    e    ? 

We  require  that,  in  response  to  an  error  such  that 

e  +    e    ,  the  system  tend  toward   e  =  z    .   That  is,  the  state 

of  the  plant,   x^  ,  should  tend  toward   x   +  e  (static 

P  m  —    v 

stability)  . 


t 


FIGURE  6.   PLANT  TRAJECTORIES  FOR  DECREASING  THE  ERROR 
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If   e  >  e    ,  we  require   e  <  0  : 


e  =  (x   -  x  )  sgn  (x   -  x  )  <  0    if   e>£       (17) 
*■  m    p   6    m    p'  v  J 


(Again,  the  Signum  function  has  no  derivative  since  the 

plant  state  trajectory  does  not  cross  that  of  the  model.) 

We  have   x  *  =  x   (6)   so, 
p     m        ' 


e  =  (x  *  -  x  )  sgn  (x   -x)<0    if   e  >  e        (18) 
v  p     pJ      °      y    m    p' 


The  two  cases  are   x   >  x    and   x   <  x   : 

m    p        m    p 


x>x   =>x*<x    \ 

m     p       P      P   j 

\      for   e  <  £  (19) 

x     x   =>  x  A  ^  x     J 

m  <   p       p   >   p   7 


These  two  cases  are  shown  in  Figure  6. 
From  equation  (1) : 


x   =  £-  (x   ,  t)  +  g.  KA  e 
p     1    p        °1 


(20) 


V  "  £1  <-*?*    ,  t)  ♦  gl  K*  e 


In  determining  whether   x  *  <  x    or   x  *  >  x   ,  we 
know   e  >  e  ,  and  we  can  examine   f,  (x   ,  t)   to  see  if 
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it  increases  or  decreases  with   x   .   But  the  sign  and 
magnitude  of   K*   will  be  determined  by  the  model  being 
followed  (equation  (15)). 

From  equation  (20)  we  have 


x   =  x  *  = 


=  [*!  Up  >  t)  -  f2  (x  *  ,  t)]  +  gl  K*  (e  -  e)  (21) 


From  equation  (19) ,  with   e  >  e  , 

First  Case:   x   >x  =>  x  *  <  x  => 

m p     p     p 


[*!  Cx   ,  t)  -  f1    (x  *  ,  t)]  +  §1  K*  (e  -  e)  >  0 


or   [f1  (xp  ,  t)  -  f1  (xp*  ,  t)]  >  -  g±   K*  (e  -  e)      (22a) 


Likewise , 

Second  Case:   X   <  x  => 
m    p 


[fx  (xp  ,  t)  -  f1    (xp*  ,  t)]  <  -  §1  K*  (e  -  e)      (22b) 

For  any  given   x   ,  x   .  e   and   g,  ,  all  the  quantities 
/  &        p  '   m  "1  '  n 

in  equations  (22a)  and  (22b)  except   K*   are  determined. 
Since  we  wish  our  descriptions  of  the  plant  and  the  model 
to  be  arbitrary,  we  cannot  assure  that  equations  (22a)  and 
(22b)  will  hold  for  all  cases. 
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This  is  illustrated  by  Example  1,  where 


x   =  a  x   +  u 

P      P    F 


x   =  b  x   +  u 
m      m    s 


Using  equations  (22a)  and  (22b) 


x   >  x    requires 
m    p     n 


ax   -  a  x  *  >  -  K*  (e  -  e) 
P      P         v  J 


a  (x   -  x  *)  >  -  K*  (e  -  e) 
P    P 


For   x   >  x   ,  (x   -  x  *)  =  -  (e  -  e) 

m    p     p    p  J  K  ' 


since   (e  -  z)    >    0  ,  we  require 

-  a  >  -  K*   or   a  <  K* 


Similarly,  for   x   <  x   we  require 
'  '       m    p       n 


a  <  -  K* 
Thus,  if   a  <  0  ,  we  require.  |a|  <  |K*|  ,  or 
(b-a)x   +  a  £  sgn  (x   -  x  *) 


a   < 


m m p_ 


using  equation  (15).   This  places  an  unacceptable  restric 
tion  on  the  choice  of  model  and  z    .       If   a  >  0  ,  one  or 
the  other  of  the  inequalities  is  not  satisfied. 
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Using  the  values  given  in  Example  1 ,  the  value  of   K- 
at   t  =  1.5   seconds  was  found  (from  the  simulation)  to  be 
-  12.2  .   Here,   a  =  0.5  ,   a   is  not  less  than   K*   and 

we  can  expect  divergence  if   x   >  x 

r         °  m    p 

A  negative  step  disturbance  of  magnitude  0.3  was  super- 
imposed on   x    at   t  =  1.5   seconds  in  the  simulation.   The 
resulting  divergence  is  shown  in  Figure  7. 


XDOT  (PLANT)  =    0.5*X  (PLANT)  +U  (PLRNT) 
XDOT  (MODEL)  =-l.xX(M3DEU+U(ST5TEM) 
PLANT    I.C. =0.05.     MODEL     I.C.=0.0. 
EPSILON    =    0.10.       NEC    DISTURB     (-0.3) 


AT    T    =    1.50 


LEGEND 
XM  □ 

XP  © 


FIGURE  7.   DIVERGENCE  RESULTING  FROM   K*e  CONTROL 
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C.   STABILIZING  THE  SYSTEM 

Neutral  static  stability   (e  =  0)   can  be  established  by 
taking   e  =  e.    That  is,  if  an  error   (e)   is  present  in 
the  system,  we  take  that  error  as  our  new   e  .   If  the 
error  changes,  we  change   £  .   Equation  (15)  becomes 

*   -  f1M  (x-  ,  e  ,  t)    *   -  f0  (x   ,  t) 

v      _     m  IN  *•  m  m    2  v  m  n- 

°l  e  g2 

where  the  subscript   N   denotes  a  neutrally  stable  system, 
and 

f1M    (x      ,    e    ,    t)    =    f,     (x   M    ,    t) 
IN    v    m    '         '      J  1    v    pN    '      } 


x  XT   =   x      -    e    sgn    (x      -    x  ) 

pN  m  °      v    m  p^ 


Since      uT,xr   =    Kx.   e 

KiN  IN 


K    -    flN    txm    '    e    '    t}         *m    -    f2    (xm    '    tJ  .... 

uo  =  i g, (24) 


At  this  point  we  are  no  longer  computing  a  gain,  but 

directly  synthesizing  a  control  modification  using 

x   ,  x   ,  e   and  the  descriptions  of  plant  and  model 
m  '   m  '  *  * 

dynamics . 

Equation  (24)  gives  the  control  which,  when  added  to 
the  system  control,  will  cause  the  plant  to  follow  the 
model  with  constant  error.   If  the  error  present  in  the 
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system  should  change,  the  plant  will  follow  the  model 
with  the  new  error  held  constant.   This  is  illustrated  by 
way  of  Example  1 .   Equation  (16)  becomes 


u™  =  (b  -  a)  x   +  a  e  sgn  (x   -  x  ) 

KN         *      m        &   v  m    pJ 


(25) 


XOOT  (PLONT!  =  C.|5*X  [PLANT]  +U  (PLANT) 
X DOT  (MODEL)  =-1  .«X  (MOOED  +U  .'  SYSTEM) 
P^ONi  l.C. =0.0.1  MODEL  I.C.=0.Q. 
FLOATING  ERROR  CONTROL  APPLIED 
PLANT  OUTPUT  DISPLACED 


LEGEND 
XM  □ 

XP  O 


9.00  0.30  0.31}  0.90 


1 . 20  1 . 50 

TIME 


FIGURE  8.   CONSTANT  ERROR  CONTROL  SIMULATION 
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The  simulation  was  run  with   x   being  disturbed  by 

P  . 

step  inputs  at   t  =  1.0   and   t  =  2.0   seconds  as  shown 
in  Figure  S. 

Note  that  the  error  is  constant  even  with   e  =  0  . 
This  does  not  alter  the  discussion  following  Figure  1 
regarding  steady  state  zero  error,  because  we  are  no 
longer  calculating  a  gain  as  shown  in  Figure  1.   The  present 
form  of  the  system  is  shown  in  Figure  9. 


Ui 


-£► 


Mooeu 


Xm 


— £— ,/ 


ul 


fc»*i 


COMTGOU 


^  PLAtOT 


Xp 


y? 


FIGURE  9.   FORM  OF  THE  SYSTEM  WITH  CONTROL  SYNTHESIS 

At  this  point,  if 

A.  The  physical  pant  is  perfectly  described  by  the 
plant  equation,  and 

B.  The  plant  and  the  model  have  the  same  initial 
conditions   (e(0)  =  0)  ,  and 
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C.  There  are  no  external  disturbances  (or  if  the  mean 
disturbance  is  zero) ,  then  the  control  defined  by  equation 
(24)  will  yield'  zero  error  (or  zero  mean  error)  . 

D.   ERROR  AND  ERROR  RATE  CONTROL 

If  an  error  exists  in  the  system,  it  is  desired  that 

it  be  reduced  to  zero  in  a  controlled  manner.   Figure  10 

shows  a  system  with  initial  error.   The  trajectory  of   x 

P 

with  constant  error  (equation  (24))  is  shown,  as  well  as 
a  desired  response  to  this  condition  (dashed  line) . 


a* 


to 

FIGURE  10.   CONSTANT  ERROR  AND  DESIRED  RESPONSE 
TRAJECTORIES  OF  THE  PLANT 


The  desired  response  may  be  obtained  by  observing  that 
the  plant  may  be  made  to  follow  any  model  with  constant 
error.   That  is,  if  we  can  find  another  model  (to  be  called 
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the  control  model)  which,  if  followed  with  constant  error, 
will  cause  the  plant  state  to  return  to  the  model 
trajectory,  the  desired  response  is  obtained. 

The  trajectory  of  the  control  model  is  shown  in 
Figure  11. 


*   ▲ 


CM 


-> 


FIGURE  11.   CONTROL  MODEL  TRAJECTORY 


The  notations   e    and   e     represent  the  errors 

m        cm    r 

between  the  plant  and  the  model,  and  between  the  plant  and 

the  control  model,  respectively.   Note  that   e     is  held 

cm 

constant,  and  that  the  trajectory  of   x    causes   e    to 

J      ;       cm  m 

diminish . 


From  Figure  10  the  nature  of   x     is  seen:   If   e    is 

°  cm  m 

0  ,  then   x     is  the  same  as  x    (except,  perhaps,  for 
'        cm  m       r    r     r 

initial  conditions).   If   e  is  non-zero,  then   x 

m  cm 

differs  from   x   by  an  amount  proportional  to   e 

m  f      t-  m 


A  control  model  \\rhich  achieves  the  desired  response  may 
be  written: 


xcm  =  *m  A  (V  sSn  (xm  '  V  ^^ 


where  A  (e  )  >  0   and 


A  (e  )  =  0   if 

v  m 


e   =  0 
m 


The  Signum  function  assures  that 


o 


£  m  <  x     if    x   >  x    and 

cm    m  p    m 


X- >  x     if    x   <  x   ,  as  required 

cm    m  p    m        n 


From  equation  (2) 


x    =  f-  (x   ,  t)  +  g0  u   +  A  (e  )  sgn  (  x   -  x  1   (27) 
cm     2  v  m  '  J         &2   s      v  nr  a      K      m    py   v  J 


so 


u   -  *cm  "  f2  (xm  ;  t}  -  &  (em}  SS"  (xm  '  V    ,,„, 

g2 
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The  additive  control  input  necessary  to  follow  the 
control  model  with  constant  error  is  then  (from 
equation  (24) )  : 


x    -  f,.,  (x    ,  e    ,  t) 
cm IN___cm_J   cm 


(29) 


x   -  f~  (x   ,  t)  -  A  (e  )  sen  fx   -  x  ) 

cm     2  v  m  '  J  v  m   to    m    p 


In  equation  (29)   f,KT  (x    ,  e    ,  t)   is  defined  as 

n       v      J         IN  v  cm  '   cm  '  J 

the  function  evaluated  at 


x   =x    -e    sgn  (x    -  x  ) 

p     cm    cm   &    cm    p 


This  expression  mav  be  shown  to  be  independent  of   x 

r  J  r  cm 

by  observing  the  six  possible  relationships  between 

x   ,  x    and   x    and  writing  the  expression  for   e 
p  '   m        cm  &       r  cm 

for    each: 


x>x>x  e        =x-x        +e 

p  m  cm  cm  m  cm  m 


x>x        >x  e        =   x      -    x        +e 

p  cm  m  cm  m  cm  m 


x     >   x     >   x  e        =  x     -    e        -e 

m  p  cm  cm  m  cm  m 


x>x        >x  e        =-x+x„+e 

m  cm  p  cm  m  cm  m 


x        >x>x  e        =    -    x      +   x        -e 

cm  p  m  cm  m  cm  m 


x        >   x     >    x  e        =    -    x     +   x        +e 

cm  m  p  cm  m  cm  m 
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From  which,  for  all  cases 


e    =  x    sgn  fx    -  x   ) 

cm    cm  °      K    cm    p7 


x  sgn  fx 
m   to 


cm 


xj 


+    em  sgn  (x   -  x  )  sgn  (x    -  x  ) 
m   to    m    p    &   ^  cm    pJ 


then 


e   sgn  (x    -  x  ) 

cm   °   v  era    pJ 


x^   -  x   +  e   sgn  (x   -  x  ) 

cm    m    m   °   v  m    p 


so 


'cm 


e    sgn  (x 

cm   5 


cm    p^ 


x   -e   sgn  (x   -  x  ) 
m    m   &    m    pJ 


The  expression  for   uK   is  then 


x 


UK 


cm 


flN  (xm 


m 


t) 


§1 


(30) 


cm 


-  f9  (x_  ,  t)  -  A  (O  sgn  (x 


m 


nr 


m 


x  ) 


where   f ,  ,T  (x   ,  e   ,  t)   is  the  function   f,   evaluated 
IN  v  m  '   m  1 


at   x   =  x   -  e   sgn  (x   -  x  ) 
p    m    m   &   v  m    pJ 


In  equation  (30)  ,  if   g,  =  g?      then  u„      does  not 


depend  on   x 


cm 


In  any  case,  it  is  simpler  at  this  point  to  replace 

the  last  term  in  equation  (30)  by   u„   and  calculate 

directly   u   =  uT,  +  u   ,  or 
7    p    K    s  ' 
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*cm  flN  (xm  ;  em  »  t} 
P  ~  g1  (3D 


The  selection  of  the  function  A  (e  )   appears  to  be 
arbitrary  with  unconstrained  control.  The  plant  following 
the  control  model  in  exactly  the  desired  manner  in  each 
simulation  run.   The  form  of  A  fe  )   chosen  for  most 
simulations  was   A  (e  )  =  [exp  (w  .  e  )  -  1]  ,  where   w 
was  selected  to  vary  the  speed  of  the  response. 

The  application  of  this  form  of  error  control  is 
illustrated  by  continuing  Example  1 : 


x   =  a  x   =  u 
P      P    P 


a  f   b 


x   =  b  x   +  u 

m      m    s 


From  equation  (26),   xcm  =  b  cm  +  us  +  A  (em)  s£n  (xm  "  xp) 


From  equation  (30),   u 


K 


x    a  (x   -  e   sgn  (x   -  x  )) 
cm    y   m.  m   °    m    p 


x    -  b  x   -  A  (e)    sgn  (xm  -  x^) 
cm ni mJ      &    m p_ 
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The  expression  for   u~   simplifies  to 


uv   =  b  *m   +  A  (ej  sgn  (xra  -  xj 


in 


nr 


m 


a(x   -e   sgn  (x   -  x  ) ) 


selecting   A  (e  )  =  [exp  (3  e  )  -  1]  . 

The  problem  was  simulated  with  the  values   a  =  +  0.5  , 
and   b  =  -  1.0   as  before.   The  plant,  model,  and  control 
model  had  the  same  initial  conditions  (zero).   A  step 
disturbance  was  imposed  on  the  plant  at   t  =  0.5   second. 
The  response  is  shown  in  Figure  12: 


XDOT  (PLANT)  =0.5*X  'PLANT)  +U  (PLANT) 
XDOT  (MODEL i  =-1  -  *X  (MODEL)  +U  (SYSTEM: 
DELTA  '=  EXP  (3.  *  (EM-EDS)  )  -1.0    ri_n 
POS  1.0  DISTURBANCE  AT  T=0.5  (STEP) 


■  ■  m  m — »  ■  » 


LEGEND 
XM       □ 
XP       O 


FIGURE  12.   SIMULATION  USING  CONSTANT  EPROk 
MODEL  FOLLOWING  CONTROL 
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The  system  at  this  point  is  most  simply  described 
from  equation  (31)  as  a  feedforward  control  synthesizer  with 
error  feedback,  as  shown  in  Figure  13: 


U^ 


WoT>El_ 


&u*a 


T 


*CM 


JLtt_ 


>" 


Up 


&HP\J\kJT 


■X§)-4- 


-M 


1 — ► 


FIGURE  13.   CONTROL  SYSTEM  CONFIGURATION  WITH 

CONSTANT  ERROR  MODEL  FOLLOWING  CONTROL 
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V.   SUMMARY 

The  method  of  constant  error  model  following  control 
is  summarized  as  follows: 
Given  the  plant  described  by 


X   =  f,  (x   ,  t)  +  g,  u  (1) 

p     1    p         &1   p 


and  the  model 

x   =  f0    (x   ,  t)  +  gn    u  (2) 

m     2  v  m  '   '    62   s 

Define  the  control  model 


x    =  x   +  A  (e  )  sgn  (x   -  x  )  •         (26) 
cm    m      v  nr   to   v  m    p^ 


Select      A    (e    )       such    that 


A    CO         0 
m     — 


A     (0)    =    0 

A    (e    , )    >    A    (e    -J       if      e    ,    >    e    ~ 
v   mlJ  ^    mlJ  ml  m2 


Solve    equation    (1)    for      u      : 


x      -    £-     (x      ,    t    ) 

_2 Lr__iU (32) 


P 


CT 
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In  equation  (32)  substitute 


where 


x      =   x 

p  cm 


f,     (x„    ,    t)    =    flv    (x      ,    e      ,    t) 
1        p  IN    v   m    '      m  J 


fiw    Cxm    ,    em    ,    t)    =   f,    (x„    ,    t)        with 


IN    v   m    '      m 


1    ^p 


x      =   xxr 

p  N 


e      sgn    (x      -    x   ) 


The    result 


u 


-    f,,T    fx      ,    e      ,    t) 

cm  IN    ^   m    '      m 


(51) 


Is  the  plant  input  which  makes  the  plant  states  equal  to 
those  of  the  model,  and  restores  this  condition  in  the 
presence  of  error. 


VI .   CONCLUSIONS  AND  RECOMMENDATIONS 

A.  It  is  concluded  that  constant  error  model  following 
control  achieves  the  desired  goal  of  error  control  for  the 
systems  described  in  this  paper.   It  offers  the  advantage 
of  perceptual  simplicity  in  that  the  designer  can  visualize 
exactly  the  effect  of  his  control  method  on  the  system 
response.   It  affords  flexibility  in  that  the  choice  of  the 
restoring  function  A  (e  )   is  arbitrary  within  the  few 
constraints  mentioned. 

B.  The  following  areas  of  future  research  are  suggested: 

1.  Development  of  the  theory  with  application  to  higher 
order  systems  with  non-additive  controls. 

2.  Investigation  of  the  response  of  systems  in  which 
the  physical  plant  is  not  accurately  described  by  the  plant 
equations . 

3.  Application  of  the  theory  to  practical  problems. 
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APPENDIX  A:   EXAMPLE  OF  APPLICATION  OF  THE  METHOD  TO 
A  FIRST  ORDER  NONLINEAR  SYSTEM  WITH 
TIME  VARYING  COEFFICIENTS 


1.   Assume  the  Plant  is  given  by 


x   =  tx  2  +  2 up  (Al) 

P     P      F 


and  the  model  by 


x   =  -  (t  x  )1/2  +  u  (A2) 

m         nr       s 


define 


x    =  x   +  A  (e  )  sgn  (x   -  x  )         (A3) 
cm    m      ^  nr  &      v  m    pJ  "•      J 


then 


x   -  tx 

u   =   E      P  (A4) 

p        2 


into  (A5)  substitute 


x   =  x  (A5) 

p    cm 


x   =  x   -  e   sgn  (x   -  x  )  (A6) 

p    m    m  5   v  m    p^ 
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then 


1  n 
(txm)   -  +  us  +  A(emj   sgn  (xm  -  x)    -  t  (^  -   ^  sgn  (xm  -  xpJ 

UP  '  2 


(A7) 
The  system  was  simulated  using  the  Continuous  System 
Modeling  Program  (CSiMP)  on  the  IBM- 5035.   Time  histories 
of  the  response  following  a  disturbance  at   t  =  0.5 
second  are  shown  in  Figures  Al  through  AS.   The  system 
control  input  is  a  unit  step  at   t  =  0   in  all  cases. 
Shown  also  are  the  plant  controls  generated  in  each  case 
The  four  cases  were  the  same  except  for  the  choice  of 

A.   A(e  )  =  5e   (Figures  Al  and  A2) 
m      m    ° 

A(e  )  =  exp  (e  )  -  1  (Figures  A5  and  A4) 
A(e  )  =  exp  (2e  )  -  1  (Figures  A5  and  A6) 
A (e  )  =  exp  (5e  )  -  1  (Figures  A7  and  A8) 
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XOOT  (PLANT)  =  TxXP*x2  +  2xl!P 
-  XOOT  (MODEL)  =-SQRT  (T*XM)  +US 

PLANT  I.C. =0.0.  MODEL  I.C.=0.0, 
DELTA  IS  5.0xEM 


LEGEND 
XM  □ 

XP  O 


FIGURE    Al. 


MODEL   AND    PLANT    RESPONSES, 
A(eJ    =    5e 


nr 


m 
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XOOT  (PLPNT)  =  T*XP**2  +  2*UP 
„  XDOT  (MODEL)  =-SQRT  (T*XM)  +US 
*,  PLANT  I.C. =0.0,  MODEL  I.C.=0.0, 

DELTA  IS  5.0xEM 


LEGEND 


UP 


FIGURE  A2  :   PLANT  CONTROL,  A(e  )  =  5e 

m      m 
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XQOT  (PLRNT)  =    TxXP*h2  +  2hIIP 
-    XQOT  (MODEL)  =-SCRT  'TkXM) ^U5 
*,  PLANT    I.C. =0.0,     MODEL     I.C.=0.0. 

DELTA    IS    EXP  (Eft)  -1 


LEGEND 
XM  □ 

xp  o 


FIGURE    A3 


MODEL  AND  PLANT  RESPONSES,  A ( e  ) 


exp(em) 
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XDOT(PLGNT)  =  TkXPx*2+2*UP 
XCOT  (MODEL)  =-SQRT  !T*XM)  +US 
PLANT  I.C.=Q.C.  MODEL  I.C.=0.0. 
DELTfl  IS  EXP  (EM)  -1 


LEGENO 


UP 


FIGURE  A4:   PLANT  CONTROL 


A(0 

m 


exp 


(■») 
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XDCT  (PLANT)  =  T*XP**2+2*UP 
-  XDQT  (MODEL)  =-SQRT  (T*XM)  +  US 
a  PLfiNT  I.C. =0.0.  MODEL  I . C. =0 

DELTR    IS     (EXP  (2.0*EM)  -1.0) 


«1 


LEGEND 
XH  □ 

XP  O 


=3.00 


1.80  2.00 

TIME 


FIGURE    A5:      MODEL    AND   PLANT    RESPONSES 


^   m 


=    exp(2 


H? 


-    1 
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XDQT  (PLANT)  =  T*XPkk2+2kUP 
_  XQOT  (MODEL)  =-SQRT  (TkXM)  +U5 
s  PLRNT  I.C. =0.0.  MODEL  I  . C. =0. 
IS  (EXP  (2.0*EM)  -1.0) 


T"B€£TR 


LEGEND 


UP 


o.oo  o.io  o.ao  1.20 


1.80  2.00 

TIME 


2.U)  2.80  3.20  3.80  U.00 


FIGURE    A6:       PLANT    CONTROL,    A(e   )    =    exp    (2    e  )     -    1 
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XDOT  (PLANT)  =  T*XP**2+2xUP 
_  XDOT  (MODEL)  =-SQRT  (T*XM)  +US 
K   PLANT  I.C. =0.0.  MODEL  I.E.: 

DELTA  IS  EXP  (3*EM)  -1 


0.0, 


LEGENO 
XM       □ 
XP       O 


FIGURE  A7:   MODEL  AND  PLANT  RESPONSES 


AOJ 
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XD2T CPLPNT) =    T*XP**2+2*UP 
XD3T  (MODEL)  =-S3RT  [T*XM1  +US 

PLANT  I.C. =0.2.  MODEL  I.C.=0.0. 
DELTA  IS  EXP  (3*EM;  -1 


LEGENO 


3.30  3.U  3.90 


L.SO  2.30  2.*0  2.30 
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FIGURE    A8 
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exp    (5    em) 


49 


LIST  OF  REFERENCES 


Motyka,  Paul  R. ,  "Variable  Stability  Simulation 
Techniques  for  Nonlinear,  Rate  Dependent  Systems," 
Proceedings  of  the  1972  Joint  Automatic  Control 
Conference . 

Rynaski,  E.  G.,  "Adaptive  iMultivariable  Model  Following," 
Proceedings  of  the  1980  Joint  Automatic  Control 
Conference . 

Rynaski,  E.  G.,  etal,  Preliminary  Design  and  Training 
Application  Analysis  of  a  YT-2B  In-Flight  Simulator, 
report  prepared  by  the  CALSPAN  Corporation,  Buffalo, 
New  York  for  the  Naval  Air  Development  Center, 
Warminster,  Pennsylvania,  under  contract  N62269- 73-C- 0937 , 
January  1974  (CALSPAN  Report  Number  AK- 5362  - F- 1)  . 


50 


BIBLIOGRAPHY 


Motyka,  Paul  R.  ,  Rynaski,  Edmund  G.,  and  Reynolds,  Philip  A 
Theory  and  Flight  Verification  of  the  TIFS  Model- 
Following  System,  Journal  of  Aircraft,  V.9,  No.  5, 
pages  347-353,  May  1972. 

Landau,  Joan  D. ,  Adaptive  Control  ,  The  Model  Reference 
Approach,  Marcel  Dekker,  Inc . ,  1979. 


51 


INITIAL  DISTRIBUTION  LIST 


No.  Copies 


1.  Defense  Technical  Information  Center  2 
Cameron  Station 

Alexandria,  Virginia   22314 

2.  Library,  Code  0142  2 
Naval  Postgraduate  School 

Monterey,  California   93943 

3.  Department  Chairman,  Code  6  7  1 
Department  of  Aeronautics 

Naval  Postgraduate  School 
Monterey,  California   93943 

4.  Dr.  Marie  D.  Hewett  1 
HR  Textron,  Inc. 

24  8  5  McCabe  Way 

Irvine,  California   92714 

5.  Mr.  E.  G.  Rynaski  1 
Arvin/Calspan  Advanced  Technology  Center 

P.  0.  Box  40  0 

Buffalo,  New  York   14225 

6.  Commander  W.  C.  Durham,  USN  2 
General  Delivery 

Dahlgren,  Virginia   22448 


52 


Thesis 

D8U 

c.l 


Durham 

Error  control  in 
model  following  control 
systems  using  constant 
error  model  following. 


2C6 


\  ji    *\  r\ 


ob 


Thesis 
c.l 


Durham 

Error  control  in 
model  following  control 
systems  using  constant 
error  model  following. 


